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Abstract 

Two-dimensional (2D) periodic electric modulations of a 2D electron gas split each Landau level 
into the well-known butterfly-type spectrum described by a Harper-type equation multiplied by an 
envelope function. This equation is slightly modified for 2D magnetic modulations but the spectrum 
remains qualitatively the same. The same holds if both types of modulations are present. The 
modulation strengths do not affect the structure of the butter fly- type spectrum, they only change 
its scale or its envelope. The latter is described by the ratio a of the flux quantum h/e to the 
flux per unit cell. Exact numerical and approximate analytical results are presented for the energy 
spectrum as a function of the magnetic field. For integer a the internal structure collapses into a 
band for all cases. The bandwidth at the Fermi energy depends on the modulation strength, the 
electron density, and, when both modulations are present, on the phase difference between them. 
In the latter case if the modulations have a tt/2 phase difference, the bandwidth at the Fermi 
energy is nearly independent of the magnetic field and the commensurability oscillations of the 
diffusive contribution to the resistivity disappear. 

PACS numbers: 73.20.At; 73.61.-r; 73.43. Qt 
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I. INTRODUCTION 



In the last fifteen years the magnetotransport of the two-dimensional electron gas (2DEG), 
subjected to periodic potential modulations, has attracted considerable experimental and 
theoretical attention. For one- dimensional (ID) modulations novel oscillations of the magne- 
toresistivity tensor have been observed at low magnetic fields B, distinctly different in 
period and temperature dependence from the usual Shubnikov-de Haas (SdH) ones observed 
at higher B. These oscillations reflect the commensurability between two length scales |2j: 
the cyclotron diameter at the Fermi level 2R C = 2^2irn e £ 2 , where n e is the electron density, 
I the magnetic length, and a the period of the potential modulation. 

The situation is similar but less clearcut for 2D electric modulations from both a the- 
oretical 3], 0] and an experimental ^] point of view. In general, for 2D modulations a 
tight-binding treatment shows that each Landau level exhibits the well-known butterfly- 
type spectrum described by a Harper-type equation. For sinusoidal modulations the energy 
spectrum resulting from the numerical solution of this equation shows, when the energy E 
is measured in units of V^Fniu^), n = x, y, with the modulation potential and F n {u^) the 
Laguerre polynomial, a nontrivial structure as a function of a = $/$o ; where $0 = h/e is 
the flux quantum and $ the flux per unit cell: for a = p/q rational, with q,p integers and 
relative prime, each Landau level is split into p subbands Q|. To our knowledge, despite 
various efforts []| there is no conclusive experimental evidence for this structure yet presum- 
ably because the small gaps between these subbands are closed due to disorder in samples 
of not exceptionally high mobilities. By neglecting those small gaps we recently accounted 
for the observed [^] oscillations in the amplitude of the commensurability oscillations, as 
a function of the flux through the unit cell of the 2D modulation lattice, that could not be 

n 

fully explained by earlier semiclassical theories [3|. The high resistivity peaks observed for a 
integer were due to collisional contributions and show up only for short-period modulations. 

In view of recent work on ID magnetic modulations it is of interest to present a 

tight-binding treatment of 2D magnetic modulations in order to see the similarity and the 
differences between the former and the latter as well as between 2D electric and magnetic 
modulations. So far we are aware only of the work of Ref. 12 which considered only a 2D 
magnetic modulation and didn't show any spectra as a function of a = q/p. In addition, 
we would like to present such a treatment when both types of modulations are present. 
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The reason for considering the latter case is that from experimen tally known methods of 
producing a magnetic modulation one expects that the magnetic [10( or superconducting 
stripes periodically placed on top of a 2DEG, also act like electrical gates and induce 
an electric modulation of the 2DEG. As shown in our previous work on transport 13] and 
the related band structure [1J], the phase difference between the two ID modulations can 
have important consequences and even suppress the Weiss oscillations; this prediction was 
later confirmed experimentally ^(|. We show here that similar effects must occur for 2D 
modulations with a tt/2 phase difference since the bandwidth at the Fermi energy is nearly 
independent of the magnetic field and the commensurability oscillations of the diffusive 
contribution to the resistivity disappear. Further, as for ID modulations 13], the bandwidth 
at the Fermi energy for a purely electric modulation oscillates, as a function of the magnetic 
field, in antiphase with that for a purely magnetic one. In addition, we show how a small 
broadening can close the very small gaps in the density of states for a = q/p rational. 

In the next section we derive the relevant Harper-type equations and solve them exactly, 
for a rational, and analytically, for a integer, to obtain the one-electron eigenfunctions and 
eigenvalues. We also present the density of states and the bandwidth at the Fermi energy. 
Various numerical results are presented in Sec. Ill and concluding remarks in Sec. IV. 



II. FORMALISM 



We consider a 2DEG in the (x, y) plane under a perpendicular magnetic field B, modu- 
lated weakly and periodically along the x and y directions, and accompanied by a 2D electric 
modulation of the same period. In Fourier space the total electric and magnetic fields can 
be expressed, respectively, by a vector potential 

A(r) = x A l^' r + y[Bx + J2 A l el&T ) W 
g V g / 

and by a scalar potential 

#r) = £> g e**\ (2) 
g 

Here a x = 2ir/K x and a y = 2ir/K y are the lattice periods of the imposed modulations and 
g = x.n x K x + yn y K y is the corresponding reciprocal lattice vector with n x , n y integers. For 
V • A = 0, which can be satisfied in most practical situations, we have A x g x + A y g y = and 



3 



the magnetic field reads 

B=Vx A = z(5 + ^5|e ig ' r ), (3) 

g 

with Bg = ig x A^ — ig y A x % the Fourier component of the magnetic modulation. 

Neglecting the spin of the electron and the Zeeman term, the one-electron Hamiltonian 
of the 2DEG reads 

H = [p + eA(r)] 2 /2m* + V(r) 

= p 2 x /2m* + (p y + eBx) 2 /2m* + ^ e 4g r (^p, + A\ Vy + eBxA y s ) 

g 

+2^E^ eig ' r +E^ eig ' r ' ( 4 ) 
g g 

where p is the momentum operator, D s = S g '(^ g - g '^g' + ^g-g'^g') * s ^ ne sc L uare term of 
the vector potential, and V g = —e<f) s is the Fourier transform of the electric potential. 

In the absence of both modulations, i.e., for B z % = and g = 0, the system has an energy 
spectrum of Landau levels E n = hu c (n + 1/2). In the Landau gauge, the corresponding wave 
functions are 

(r\n, k y ) = ip Ujky = L~ 1/2 exp(zA; y y)</> n (x + x ) (5) 

with xq = l 2 k y , I 2 = h/eB, and <ft n (x) the normalized harmonic oscillator function. Con- 
straining the center of the cyclotron orbit — Xq to be within the sample and using periodic 
boundary condition along y, gives the degeneracy of the Landau levels as L x L y /2irl 2 = N, 
with L x and L y the sample dimensions. We choose N to be an integer. 



A. Tight-binding treatment 

In the presence of weak ID modulation, when coupling between different Landau levels due 
to modulation is negligible, the Landau levels are broadened into bands and their bandwidth 
oscillates with the uniform field B. The eigenstates of the system and its energy spectrum 
can be calculated using perturbation theory or following a tight-binding treatment ^| . 

For 2D modulations, the situation is much more complicated and in general the eigen- 
states cannot be the same as those without modulation. A tight-binding method |3] used 
by Hofstadter [f| to treat a 2DEG with strong electric modulation but weak perpendicular 
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magnetic field results in a "butterfly" -type energy spectrum. Neglecting coupling between 
different Landau levels, we can express the new eigenstates as linear combinations of the 
states without modulation for each Landau level [17[. In the following, we will show that 
using the same method the electron states in the presence of both electric and magnetic 
modulations can also be obtained. Notice, however, that in this method the parameter a is 
replaced by its inverse, [18j i.e., in what follows we will have a = $o/^ — l/p- 

For weak modulations, i.e., for B^<^ B and V g <C fru c , we neglect the coupling between 
Landau levels. Then from Eq. (4) we obtain, for the nth Landau level, the matrix elements 

(n, k' y \H\ n, ky) = EJ k ,, ky + £ U^^+^H^ , ky+9y} (6) 

g 

where U g is the effective potential 

U e = V g F n (Q) + (eh/m*)B z s G n (Q) + (e 2 /2m*)D g F n (Q), (7) 

with F n (x) = exp(-x/2)L n (x), G n {x) = -(d/dx)F n (x), Q = l 2 (g 2 x + £f 2 )/2, and L n (x) the 
Laguerre polynomial of order n. 

^From Eq. © we see that if only the modulation along the x direction is present, i.e. 
for 9v = 0, there r esu. tS on.y a .evel broadening for each state but no coup.ing between 
different states. As noticed earlier [9|, [15|], the modulation along the y direction, however, 
introduces a coupling between the states (n, k y \ and (n, k y + K y \. This leads to the folding 
of the energy band and the formation of a Brillouin zone along the y direction. Then we can 
assume a wave function of the form: ||n, () = J2\ c a(^, C) l n > ky + ^Ky) with A = 1, 2, ...M 
for integer M = Na y /L y . The secular equation H \\n, () = £ n ^ \\n, () then becomes 

£ {H\/\ - (£ nX - E n )6 yx } c A (n, C) = 0, (8) 

A 

with Hy\ = J2 g Sy t x+ gy /K y U g exp[—il 2 g x (2k y + (A' + X)K y ) /2]. For sinusoidal modulations 
and rational values of a Eq. (8) becomes an extension of Harper's equation, see below. 
For a rational value of 

a = $ /$ = 2nl 2 /a x a y = q/p (9) 

with q and p integers and relative prime, the M-dimensional matrix H\ t \> satisfies the 
periodic condition Hx+ Pt \>+ P = H X y. Due to the peculiar character of the harmonic function 
centered at x = l 2 k y , this corresponds to an effective spatial periodicity of the Hamilitonian 



with period pl 2 K y = qa x and Bloch-type extensive states along the x direction emerge as 
a general property of the superlattice. It is convenient [l2| to introduce a new Bloch-type 
basis \\\n,s,k x ,k y ) with a new quantum number k x , the momentum of electrons along the 
x direction, 

N x 

|| \n, s, k x , k y ) = (Ay 1 / 2 e~ a2k ^y^ p+s ^ \n, k y + {(3p + s)K y ) , (10) 

13=1 

where N x = L x /qa x is assumed an integer by choosing L x . 
In this basis the eigenstate of the Hamiltonian reads 

p 

I \ n -> ji kxi ky)) = u s (n, j, k x , ky) || \n, s, k x , ky) . (H) 

s=l 

Since u s {...) obeys the relation u s+p (n,j, k x , k y ) = u s (n,j, k x , k y ), it is a periodic function of 
s and the M x M Hamiltonian matrix is reduced into a p x p one with elements 

[h { v\k)] s , tS = ^2slf y/Ky] s ,_ s U ntS exp[-il 2 (g x ky - g y k x + sg x K y + g x g y /2], (12) 
g 

and 

{p) J 1 if [g y /K y ] = s' - s or s' - s + p 

%y/Ky],S'-S ~ \ . ' y 16 ) 

I otherwise 

[g y /Ky] is g y /K y modulo p, k x e [0,2ii/qa x }, k y G [0,2Tr/a y ], and s = 1,2, ••• ,p. Each 
Landau level is split into p minibands with energy denoted by S^l(k x , k y ). 

In this paper we will discuss in detail the energy spectrum and the density of states 
(DOS) of a 2DEG in a constant magnetic field B with a weak magnetic modulation B = 
[em*uj x cos(K x x) + em*uj y cos(K y y)]z and a weak electric modulation of the potential V = 
V x cos(K x x + ip x ) + V y cos(K y y + {p y ) with <p x , ip y the phase differences between the respective 
components of the potential. The corresponding effective potential is found to be 

U n = 22 U n , s e ls ' r = V x F n (u x ) cos(K x x + ip x ) + V y F n (u y ) cos(K y y + ip y ) 
g 

+Hu x G n (u x ) cos(K x x) + tvui y G n {u y ) cos(K y x) 

+e x [l - F n (4u x ) cos(2K x x)} + e y [l - F n {Au y ) cos{2K y y)}, (14) 
with = l 2 K 2 j2 and e M = m*u 2 / AK 2 ji — x,y. With 8n) s {k x , k y ) = S n ^ the corresponding 
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secular equation, Eq. (jHJ), becomes 

-^F n (4u y )[C x _ 2 + C x+2 ] + {e, + e y - e x F n {4u x ) cos[2l 2 K x {k y + XK y )\ 

+hu x G n (u x ) cos[l 2 K x (k y + \Ky)\ + V x F n {u x ) cos[l 2 K x (k y + \K y ) - ip x ] - (£ n>c - E n )} C x 

+^[tkJ y G n (uy) + V y F n (u y )e iip y}C x ^ + ^{huo y G n {u y ) + V y F n {u y )e-^]C x+1 = 0. (15) 

For a = 2l 2 K x K y = q/p rational, with p, q integers and relative prime, the exact solution 
of this difference equation is obtained numerically by rewriting it in the form of a matrix 
0], of dimension p, and solving the resulting eigenvalue problem with the details given 
above. It will be shown below. 

For a pure electric modulation Eq. (15) holds with the terms oc and oc set equal 
to zero while for a pure magnetic modulation the terms oc vanish; in either case we set 
M = 0. Comparing these two cases, we see that Eq. (15) should lead to the same overall 
spectrum if hjj^G n {u^) = V^F n {u^) provided we neglect the very small quadratic terms oc e^. 
Examples of these spectra will be shown in Sec. III. 



B. Results for a integer 

For a = 2l 2 K x K y integer the dependence of the quantity {...} in Eq. (15) on A drops out 
since A is an integer. Then Eq. (15) admits the exponential solutions C x = Coexp(z£A); the 
resulting energy eigenvalue is 

£^{k x ,k y ) = E n + e x [l- F n {4u x )cos2 V ]+e y [l- F n {4u y )cos2£} 
+Hu> x G n (u x ) cost] + t\w y G n {u y ) cos£ 

+V x F n (u x ) cos{t] - (p x ) + V y F n (u y ) cos(£ + (p y ), (16) 

where £ = k x £ 2 Kj, and r] = k y £ 2 K x . This is a simple sinusoidal band as when only one modu- 
lation is present |8| . The corresponding eigenvectors are labeled with the additional quantum 
number £ (0 < £ < 2n): \ ip n k y ^ >= CoJ2 P ex P(^^)\ n ^y + ^y >• The orthonormality 
condition gives £ = k x £ 2 K y and Cq = £{K y j L x ) x l 2 by normalization. 

If we set V x = V y = in Eq. (16), we see that the only difference from the case of a pure 
electric 2D modulation is the presence of the terms e M , which are proportional to uj 2 , on the 
first line of Eq. (16); however, these terms are very weak compared to those on the 2nd 
and 3rd line oc u^. Further, it is instructive to contrast the spectrum, given by Eq. (16), 
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for a pure electric 2D modulation (hu^ = 0) with that for a pure magnetic 2D modulation 
(Vp = 0) as done for Eq. (15). We then see explicitly that the two spectra are the same 
provided hu^G n (u^) = V^F n {u^). This means that the results of a pure electric modulation 
can be mimicked by those of a pure magnetic one and vice versa. Notice, however, that 
this equivalence does not hold for those values of for which the factors G n (u^) or F n (u^) 
vanish since they are not the same in both cases. The same conclusion was reached in Ref. 
12 with the corresponding spectrum evaluated by first-order perturbation theory. 

In the absence of modulations we have the discrete Landau levels E n . In their presence 
Eq. (16) shows that these levels are broadened into the bands £n (k x ,k y ) with explicit 
sinusoidal dependence on k x and k y . If we neglect the very small terms oc e M on the first 
line and take <p x = (p y = for simplicity, their bandwidth is given by AE n (u x ,u y ) = 
2[V x F n (u x ) + V y F n {uy) + hjj x G n {u x ) + ftw y G n (u y )} and it oscillates with the field B. The 
related velocity components [v^ = (l/h)(d£n\k y ,k x )/dkff), fi = x,y) resulting from Eq. 
(16) are 

v x = -{l 2 K y /H) [V y F n (u y ) sin(£ + <f y ) + fuj y G n (u y ) sin £ - 2e y F n (4u y ) sin2£] , (17) 

v y = -{l 2 K x /K) [V x F n (u x ) sin(?7 - ip x ) + tko x G n (u x ) sinr? - 2e x F n {Au x ) sm2rf\ . (18) 

The broadening of the Landau levels into bands that gives rise to these velocity components 
has important consequences for transport especially when the fine structure of the exact 
energy spectrum is not resolved due to disorder [8] ■ The main consequence is a non vanishing 
diffusive contribution to the conductivity which is absent when the modulations are not 
present. 

In the following we will show that we can use the expression of the energy spectrum, Eq. 
(16), and the velocity operators, Eqs. (17) and (18), for a integer to approximately estimate 
the physical properties of the 2DEG system for any value of a. 

C. The density of states 

The energy spectrum given by Sn}(k X} k y ) is qualitatively different from the unmodulated 
spectrum, given by E n , and from the corresponding ID modulation spectrum given by 
E n + F n {u x ) cos K x xq. These differences are also reflected in the density of states (DOS) 
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defined by 

D(E) = 2 S(E-S^(k x ,k y )). (19) 

For a Lorenzian broadening of width T, the 5 function in the above equation should be 
replaced by (l/7r)r/[r 2 + (E — Si P l(k x , k y )) 2 }. 

The result for a integer is obtained by replacing £n}(k x , k y ) by Sn}{k x , k y ) in Eq. (19). 
For a 2D modulation with rectangular symmetry, corresponding to Eq. (16), the DOS 
becomes 

D{E)= D J2 d£{[V x F n (u x ) + hw x G n {u x )f (20) 

n=0 J ° 

— [E — E n — {V y F n {u y ) - hw y G n {u y )} cos£] 2 }- 1/2 , 
where Dq = L y L x /7i 3 £ 2 . The radicant in Eq. (20) must be positive. 

III. RESULTS AND DISCUSSION 
A. In-phase, antiphase modulations 

The results presented so far are valid for rectangular modulations. However, as the various 
expressions and results become simpler for square modulations, we will consider only the 
latter. For a square modulation where K x = K y = K, ip x = ip y = <p, V x = V y = V, u x = u y = 
u>, the secular equation (jSJ), without the square terms oc e M due to the magnetic modulation, 
is Harper's equation with (p = if the modulations are in phase and </? = 7r if they are in 
antiphase; in these cases Eq. (8) takes the form 

+ C x+1 /2 + C x cos(A/ 2 K 2 + Kx ) - C x e n , c = 0, (21) 

where e n ^ = (£ n ,c ~ E n )/[VF n (u) ± hwG n (u)]. The resulting energy spectrum is the "but- 
terfly" spectrum en? s (k x , k y ) G [—2, 2] modulated by an envelope function determined by the 
modulations 

K) =E n + [VF n {u) ± hwG n (u)}e^(k x , ky), (22) 

where the upper (lower) sign is for in-phase (antiphase) modulations. 

The spectrum for a integer and the corresponding wave function were used previously 
as an approximation, for all values of the magnetic field, and successfully described 
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magnetotransport of a 2DEG in the presence of a 2D electric modulation in the case that 
the fine structure of the exact spectrum is not resolved. This approximation implies the 
replacement of en) s (k x , k y ) by e^\(k x , k y ) = cos£ + cos 77 in Eq. (J22J), which becomes 

k y ) ~E n + [VF n (u) ± hwG n (u)}(cos£ + cost]). (23) 

In the following the resulting spectrum will be referred to as approximate band structure or 
approximate spectrum. 

The edge of the approximate band is the envelope of the exact energy spectrum. Employ- 
ing the asymptotic expressions of the Laguerre polynomials for n>l, the corresponding 
bandwidth is estimated as 

AE n (u) ~ A(n 2 nu) ~ 1/4 7(1 + 6 2 ) 1/2 |sin(2Vmi - tt/4 + (3)\ , (24) 

with 5 = (Huj /V)y/n/u, V = V ± Huj/Au and (3 = arctan [±(V/Hu;)y/u/n], so that the flat 
band condition is j3 ~ mir + 7r/4 — 2 v /mx where n is the Landau index and u = K 2 l 2 /2. For 
1/Au <C yn/u or a > 16n7r, which is valid for typical parameters discussed in this paper, 
we can further approximate V with V and 5 with 5 = (Hon /V)^n/u. Because (3 — 7: for a 
pure electric modulation and (3 = for a pure magnetic modulation, the two modulations 
are phase shifted and their interplay will influence the results of transport experiments. 

In Fig. 1 we plot the envelope function of the energy spectrum for Landau levels n = 
to n = 10 as a function of the cyclotron energy Hlj c . For the same modulation strength 
V = Hujq = O.lmeV the ratio between the band broadenings is determined by the value 
of S, cf. Eq. (24). For a period a = 200 nm, cf. Figs. 1(a) and (b), the curve 6=1 
is located near E w 0.8meV and the broadening due to magnetic modulation is generally 
larger than that due to the electric one. Decreasing the period a, the 5 = 1 curve is shifted 
to a higher energy, higher than 3meV for a = 80nm. This is the reason why the magnetic 
band broadening is narrower than the electric one for energies lower than 3meV as shown 
in Figs. 1(c) and (d). 

Figure 2 illustrates the envelope function as in Fig. 1 with both modulations having the 
same strength V = Tiwq = 0.1 meV and being in phase. In such a case, the modulation 
of longer period results in a broader bandwidth. This is opposite to that obtained from 
the band broadening in a common superlattice and results from the presence of the strong 
perpendicular magnetic field. In a common superlattice the miniband appears as a result of 
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h(O c (meV) 



FIG. 1: Approximate band structure, for n = 0, ...10, as a function of Hlu c for square electric ((a) 
and (c)) or magnetic ((b) and (d)) modulations of long period (a = 200nm), in (a) and (b), and 
short period (a = 80nm) in (c) and (d). The crossings of the dotted curves with the Landau levels 
give the asymptotic flat-band positions and the dashed line, in panels (a) and (b), is the 5=1 
curve. 

electron tunneling between adjacent quantum wells. The miniband width decreases when 
the superlattice period increases since the tunnel coupling weakens if the other parameters 
remain the same. However, here the level broadening stems mainly from the perturbative 
correction to the Landau energy of each state by the modulation. In Fig. 3 we show the 
resulting bandwidth of the 10th Landau level as a function of the modulation period for fixed 
magnetic field B = 0.5 T and the strength parameters shown in the caption. We see that on 
the average the width increases when the period increases, say by a factor of s provided s is 
such that sa is not close to the point where the bandwidth vanishes. On the right panels of 
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FIG. 2: The same as in Fig.l for both electric and magnetic in-phase modulations of long (a) 
and short (b) period. The thin solid curves show the Fermi energy for temperature T = IK and 
electron density n e = 10 11 cm -2 . 

Fig. 3 we show the constant (£- or 77-indepedent) term of the square of the velocity, given by 
Eq. (17) or Eq. (18), vs the modulation period a for the same magnetic field. If we neglect 
the very small terms oc e M , this term is equal to [(£ 2 K y / h)V y F n (u y ) + f\WyK y G n {u y )] 2 /2. 
This term gives by far the dominant contribution to the diffusive conductivity, given by the 
standard expression cr^/(0) = (/?e 2 /f2) 5^f/c(l — fc) T (E() v ft v ii where / is the Fermi-Dirac 
function, the velocity, and r the relaxation time |s| pertinent to a state £ = {k x , k y ,n}. 
With the integrals over k x and k y giving, for this term, just a constant prefactor near the 
Fermi level, we see that this contribution is mainly determined by the bandwidth. By 
contrasting the left and right panels one can see that this on-the-average increase of the 
bandwidth with the period is mainly due to the behavior of the contribution of the magnetic 
modulation to it. This behavior of the bandwidth and the antiphase between panels (a) and 
(b) on the left, and panels (d) and (c) on the right, are directly related to those of the 
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FIG. 3: Left panels: the approximate bandwidth A of the 10th Landau level as a function of 
the modulation period a at magnetic field B = 0.5 T for modulation strengths {a)V = 0.1 meV 
and Huj = 0, (b) V = and hw = 0.1 meV, and (c) V = fouj = 0.1 meV. Right panels: the 
constant term of the square of the velocity in arbitrary units, given by Eq. (17) or Eq. (18), vs 
the modulation period a at magnetic field B = 0.5. Panels (d), (e), (g) correspond to panels (a), 
(b), (c), respectively. The two periods used in Figs. 1 and 2 are 80 nm, indicated by the triangle, 
and 200 nm. 

Laguerre polynomials and their derivatives that appear in the factors F n and G n . 

If both modulations are present, their relative contribution to the energy broadening is 
determined by the value of 5. The magnetotransport behavior of a 2DEG is directly related 
to the level broadening at the Fermi energy. Because the electric modulation dominates the 
low-energy broadening and the magnetic modulation dominates the high-energy one, the 
transition from one energy range to another can be observed if we shift the Fermi energy 
through these ranges by changing the electron density. In Fig. 4 we plot the bandwidth at 
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FIG. 4: The approximate bandwidth of the Landau level closest to the Fermi level as a function of 
the magnetic field for different electron densities n e = 0.5 (a), 1.0 (b), 2.0 (c), and 5.0 x 10 11 cm -2 
(d). The solid curves denote broadening with both modulations present and of the same strength. 
The dotted (dashed) curves show the broadening when only the electric (magnetic) modulation is 
present. 

the Fermi energy for a period a = 80nm as a function of the perpendicular magnetic field B. 
The bandwidth oscillations, resulting when only the electric modulation is present, are in 
antiphase with those resulting when only the magnetic one is present. This should show up 
in the commensurability oscillations of the magnetoresistance as it does for ID modulations 



13|. If both modulations have the same strength, the bandwidth oscillations (solid curves) 
are similar to those due to the electric one at the lower electron density (Fig. 4(a)) and 
similar to those due to the magnetic one at the higher electron density (Fig. 4(d)). The non 
smooth behavior of the bandwidth, when the field B is varied, reflects that of the Fermi level, 
as shown, for instance, in Fig. 2. Upon changing the electron density the Fermi level crosses 
different Landau levels in a non smooth manner, cf. Fig. 2. For higher densities and/or 
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FIG. 5: Exact and approximate (a integer) spectrum, in units of hu> c , vs the cyclotron energy hu c 
for the n = 2 Landau level and a period a = 80 nm. The upper panel is for an electric modulation, 
the middle one for a magnetic one, and the bottom one for both modulations present. Notice the 
continuous bands for a integer: in this case the exact, numerically obtained bandwidth coincides 
with the one obtained analytically. 

longer periods more Landau levels are occupied and the oscillations are much smoother. 

In Fig. 5 we plot the exact energy spectrum of the second Landau level, in units of fiw c , 
as a function of fru c for (a) only an electric modulation, (b) only a magnetic modulation, and 
(c) both electric and magnetic modulations all of period a = 80nm and strength V = Hlj = 
O.lmeV. The positions of a — q/p = 2,3,4,5 are marked by the dashed vertical lines. The 
exact spectrum is the same as Hofstadter's spectrum modulated by the envelope function 
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FIG. 6: Broadened DOS for the exact (thin solid curves: T = 0.01K, thick curves: T = 0.5K) and 
the a integer spectrum approximation (dotted curves: T = 0.01K, dashed curves: T = 0.5K). The 
magnetic field is such that a = 2 in (a), a = 7/3 in (b), a = 5/2 in (c). V = Hujq = O.lmeV 

(dotted curves in Fig. 5). For a rational value of a = q/p the spectrum is composed of p 
minibands which touch each other for a half-integer or p = 2. For a integer the spectrum is 
the sinusoidal band given by Eq. (23). The envelope function, given by the dotted curves, is 
the band edge of the a integer spectrum [Eq. (23)] used here for all values of a. A negligible 
shift and the modification of the spectrum, due to the terms oc e^, are not shown in the 
figures. 

In Fig. 6 we show the exact and approximate DOS for in-phase modulations and various 
values of a. The other parameters are specified in the caption. The results coincide, as they 
should, for a = 2. Notice the horizontal scale and how a small level broadening T closes the 
gaps of the exact DOS in panels (b)-(d). 



B. Out-of-phase modulations 

For arbitrary phase difference between the two modulations, the secular equation cannot 
be written simply as Harper's equation but as a more general one, cf. Eq. (15). However, the 
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energy spectrum for a integer is still the envelope of the exact spectrum and the asymptotic 
bandwidth is 

AE n (u) = 4V(7T 2 nu)- 1/4 [l + (5 2 - 1) sin 2 (2v^ - tt/4) - <5cos0cos(4v^)] 1/2 , (25) 

where 5 = (fuv/V)(n/u) 1 ^ 2 . At the Fermi energy we have 5{Ef) = akpfku / (2irV) . If 5 = 1 
and ip = 7r/2, we get a nearly constant bandwidth AV^nu)' 1 ^ 4 . 

In Fig. 7 we plot the energy of the n = 0, ...,10 Landau levels as a function of the 
magnetic field for equal-strength (0.1 meV) modulations, of period 80 nm, that are out of 
phase by tt/2. The dotted curve, i.e. 5 — 1, shows a constant bandwidth that, as explained 
above, leads to a washout of the commensurability oscillations. The thin solid curve shows 
the Fermi level at T = 5 K. 

The amplitude of the bandwidth oscillations varies with the phase difference between the 
two modulations and/or with their relative strength. In Fig. 8 we plot the bandwidth for 
equal modulation strengths, V — ftw — 0.1 meV, and different phases. As Fig. 8(b) shows, 
the amplitude can be greatly decreased when <fi = n/2 and this implies a washout of the 
commensurability oscillations. Different from the in-phase or antiphase case (Fig. 4), we 
notice that the bandwidth increases more rapidly with magnetic field and that it does not 
reach zero at specific values of B. In Fig. 9 we show again the bandwidth for modulations 
phase shifted by n/2 but with different relative strengths apart from panel (c), in which the 
two strengths are equal, and the amplitude of the oscillations minimal. 

Finally, in Fig. 10 we plot the exact and approximate spectra of the n = 2 Landau level, 
as a function of the magnetic field, for various phase differences between the two modulations 
as indicated. The approximate result, for a integer, is shown by dotted curves and is again 
the envelope function of the exact spectrum as for in-phase modulations. 

IV. CONCLUSIONS 

We studied in detail the band structure of a 2DEG in the presence of weak 2D electric 
and magnetic modulations as function of an applied perpendicular magnetic field B. The 
tight-binding description shows a Hofstadter-type spectrum with an envelope function that is 
determined by the strengths of the modulations. This envelope coincides with the spectrum 
obtained from the Harper-type equation for a = $o/^ integer, where $ is the flux per 
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FIG. 7: Band structure with both modulations present and phase-shifted by <fi = n/2. The Fermi 
level at temperature T = 5 K is also indicated by the thin solid curve. The 6 ~ 1 curve shows 
a constant bandwidth leading to a washout of the commensurability oscillations of the diffusive 
contributions to the resistivity. 

unit cell. For a integer the analytical spectrum, Eq. (16), coincides with the one obtained 
from the numerical solution of Eq. (15). As an approximation, one can use the spectrum 
for a integer for all magnetic fields if the small gaps of the exact spectrum are closed 
due to disorder [8j]. The flat-band condition and the bandwidth are then found from the 
approximate bandwidth for large n. 

For a = q/p rational, the gaps in the exact spectrum, for the usual sinusoidal 2D modula- 
tions, are small and nearly close when a small level broadening is included in the calculation 

fl 

of the DOS whether both modulations are present or only the electric one |8j. The value 
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FIG. 8: Bandwidth at the Fermi energy as a function of B for different phases (f> = tt/4 (a), <p = tt/2 
(b), and = ir (c) between the two modulations. 
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FIG. 9: The same as in Fig. 7 for different electric modulation strengths V, as indicated, and fixed 
magnetic one Hloq = O.lmeV at phase shift of ir/2. 
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FIG. 10: Exact and approximate (dotted curves) spectrum for various phase differences between 
the two modulations (p = tt/4 (a), <f> = ir/2 (b), <fi = 37r/4 (c), and 4> = ir (d). 

of T needed to close the gaps depends on the modulation strength. For instance, with 



V x = V y = 0.5 meV, a width T = 1.1 K was necessary [8|]. In both cases, T is quite small. 
Accordingly, the very interesting fine structure of the spectrum may be very difficult to 
observe. 

As discussed in the text and illustrated in Fig. 3, the spectrum for a integer and a 
pure electric modulation is in general equivalent to that of a pure magnetic modulation if 
we neglect the very small terms oc e M . The same conclusion was reached in Ref. 12 for a 
perurbative evaluation of the spectrum. The same analogy, though less explicitly, applies to 
the exact spectrum obtained from Eq. (15) and shown in Fig. 5. 

In line with our previous study ^], for other kinds of surface superlattices, e.g., hexagonal 
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or trigonal, the results are similar to those presented here even when we include cross terms 
oc V x V y cos(K x x) cos(Kyy) in the modulation potential. 

The oscillations of the bandwidth with the magnetic field B, due to only an electric 
modulation, are in antiphase with those due to only a magnetic modulation. If both mod- 
ulations are present, an additional parameter is the phase between them and the situation 
is more complex. The relative importance of the two modulations can be estimated by the 
parameter 5, which is determined by the modulation strengths and the electron density. 
The electric modulation dominates for 5 < 1 and the magnetic one for 8 > 1. By changing 
the electron density, the Fermi energy can be shifted from the 5 < 1 regime to the 5 > 1 
one. Alternatively, one can change the relative strengths between the electric and magnetic 
modulations to tune the value of 5 which subsequently will influence the magnetoresistance. 
Similar to the case of ID modulations [13], a transition between the two regimes occurs for 
5 = 1 and the commensurability oscillations of the diffusive contributions to the resistivity 
disappear if the two modulations are phase shifted by n/2. 
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